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1 Introduction 



In this talk I will discuss techniques for the computation of loop integrals, which occur in pertur- 
bative calculations in quantum field theory. Particle physics has become a field where precision 
measurements have become possible. Of course, the increase in experimental precision has to 
be matched with more accurate calculations from the theoretical side. This is the "raison d'etre" 
for loop calculations: A higher accuracy is reached by including more terms in the perturbative 
expansion. The complexity of a calculation increases obviously with the number of loops, but 
also with the number of external particles or the number of non-zero internal masses associated 
to propagators. To give an idea of the state of the art, specific quantities which are just pure 
numbers have been computed up to an impressive fourth or third order. Examples are the cal- 
culation of the 4-loop contribution to the QCD ^-function [1], the calculation of the anomalous 
magnetic moment of the electron up to three loops [2], and the calculation of the ratio of the 
total cross section for hadron production to the total cross section for the production of a /j^/j 
pair in electron-positron annihilation to order O (a^) [3]. Quantities which depend on a single 
variable are known at the best to the third order. Outstanding examples are the computation of 
the three-loop Altarelli-Parisi splitting functions [4, 5] or the calculation of the two-loop ampli- 
tudes for the most interesting 2 — > 2 processes [6-16]. For the calculation of these amplitudes, 
the knowledge of certain highly non-trivial two-loop integrals has been essential [17-19]. The 
complexity of a two-loop computation increases, if the result depends on more than one variable. 
An example for a two-loop calculation whose result depends on two variables is the computation 
of the two-loop amplitudes for e+e^ — > 3 jets [20-22]. But in general, if more than one variable 
is involved, we have to content ourselves with next-to-leading order calculations. An example 
for the state of the art is here the computation of the electro-weak corrections to the process 
e+e^ — i> 4 fermions [23,24]. 

From a mathematical point of view loop calculations reveal interesting algebraic structures. 
Multiple polylogarithms play an important role to express the results of loop calculations. The 
mathematical aspects will be discussed in this talk. Additional material related to loop calcula- 
tions can found in the reviews [25-28] and the book [29]. 

This paper is organised as follows: In the next section I review basic facts about Feynman 
integrals. Section[3]is devoted to the Mellin-Bames transformation. In section |4] algebraic struc- 
tures like shuffle algebras are introduced. Section[5]deals with multiple polylogarithms. Section[6] 
combines the various aspects and shows, how certain Feynman integrals evaluate to multiple 
polylogarithms. Finally, section |7] contains a summary. 

2 Feynman integrals 

To set the scene let us consider a scalar Feynman graph G. Fig. [T] shows an example. In this 
example there are three external lines and six internal lines. The momenta flowing in or out 
through the external lines are labelled p\, pi and p2, and can be taken as fixed vectors. They are 
constrained by momentum conservation: If all momenta are taken to flow outwards, momentum 
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Figure 1: An example of a two-loop Feynman graph with three external legs. 



conservation requires that 

P\+P2+p-i = 0. (1) 

At each vertex of a graph we have again momentum conservation: The sum of all momenta 
flowing into the vertex equals the sum of all momenta flowing out of the vertex. A graph, where 
the external momenta determine uniquely all internal momenta is called a tree graph. It can be 
shown that such a graph does not contain any closed circuit. 

In contrast, graphs which do contain one or more closed circuits are called loop graphs. If we 
have to specify besides the external momenta in addition / internal momenta in order to determine 
uniquely all internal momenta we say that the graph contains / loops. In this sense, a tree graph 
is a graph with zero loops and the graph in fig. [T] contains two loops. Let us agree that we label 
the I additional internal momenta by ki to ki. 

Feynman rules allow us to translate a Feynman graph into a mathematical formula. For a 
scalar graph we have substitute for each internal line j a propagator 

(2) 



Here, qj is the momentum flowing through line j. It is a linear combination of the external 
momenta p and the loop momenta k: 

= Qj{P:k). (3) 
nij is the mass of the particle of line j. The propagator would have a pole for p^j = mj, or 



phrased differently Ej = iyp^ + m^. When integrating over E, the integration contour has to 

be deformed to avoid these two poles. Causality dictates into which directions the contour has to 
be deformed. The pole on the negative real axis is avoided by escaping into the lower complex 
half-plane, the pole at the positive real axis is avoided by a deformation into the upper complex 
half-plane. Feynman invented the trick to add a small imaginary part id to the denominator, 
which keeps track of the directions into which the contour has to be deformed. In the following 
the i5-term is omitted in order to keep the notation compact. 

The Feynman rules tell us also to integrate for each loop over the loop momentum: 

IIl. (4) 

(271)4 
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However, there is a complication: If we proceed naively and write down for each loop an integral 
over four-dimensional Minkowski space, we end up with ill-defined integrals, since these inte- 
grals may contain ultraviolet or infrared divergences! Therefore the first step is to make these 
integrals well-defined by introducing a regulator. There are several possibilities how this can be 
done, but the method of dimensional regularisation [30-32] has almost become a standard, as the 
calculations in this regularisation scheme turn out to be the simplest. Within dimensional regular- 
isation one replaces the four-dimensional integral over the loop momentum by an D-dimensional 
integral, where D is now an additional parameter, which can be a non-integer or even a complex 
number. We consider the result of the integration as a function of D and we are interested in the 
behaviour of this function as D approaches 4. It is common practice to parameterise the deviation 
of D from 4 by 

D = 4-2£. (5) 

The divergences in loop integrals will manifest themselves in poles in 1 /e. In an /-loop integral 
ultraviolet divergences will lead to poles 1 /e' at the worst, whereas infrared divergences can lead 
to poles up to l/e^'. We will also encounter integrals, where the dimension is shifted by units of 
two. In these cases we often write 



D = 2m -2£, (6) 

where m is an integer, and we are again interested in the Laurent series in e. 

Let us now consider a generic scalar /-loop integral Iq in D = 2m — 2£ dimensions with 
n propagators, corresponding to a graph G. Let us further make a slight generalisation: For 
each internal line j the corresponding propagator in the integrand can be raised to a power \ j. 
Therefore the integral will depend also on the numbers Vi,...,v„. We define the Feynman integral 
by 

^ I f ' d^Jc " 1 

J r=i iTi^ j=i {-rj+tnjyj 



The momenta qj of the propagators are linear combinations of the external momenta and the 
loop momenta. In eq. dV]) there are some overall factors, which I inserted for convenience: ^ is 
an arbitrary mass scale and the factor /j^^ ensures that the mass dimension of eq. (|7]) is an integer. 
The factor e^^^ avoids a proliferation of Euler's constant 



Je = lim >- -Inn =0.5772156649... (8) 



in the final result. The integral measure is now d^k/iin^/^) instead of d^k/{2Ti)^, and each 
propagator is multiplied by /. The small imaginary parts /6 in the propagators are not written 
explicitly. 

How to perform the D-dimensional loop integrals ? The first step is to convert the products 
of propagators into a sum. This can be done with the Feynman parameter technique. In its full 
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generality it is also applicable to cases, where each factor in the denominator is raised to some 
power V. The formula reads: 




n 



5 1 - L X, 



i=lPi Ur(vAn \;=1 / / v^.D. \ i=i 



!=1 



Applied to eq. dV]) we have 



Y,x,P, = Y.x,{~qj + mj). (10) 

i=l i=l 

One can now use translational invariance of the D-dimensional loop integrals and shift each loop 
momentum kr to complete the square, such that the integrand depends only on kj. Then all D- 
dimensional loop integrals can be performed. As the integrals over the Feynman parameters still 
remain, this allows us to treat the D-dimensional loop integrals for Feynman parameter integrals. 
One arrives at the following Feynman parameter integral [33]: 



nr(v,) ^0 Vi=i 



;=1 



The functions U and depend on the Feynman parameters. If one expresses 

n III 
Y^Xj{~q]+m)) = -Y^Y.krMrsks+Y.'^kr-Qr-J, (12) 

j— 1 r— 1 s= 1 r= 1 

where M is a / x / matrix with scalar entries and 2 is a /-vector with fourvectors as entries, one 
obtains 

'U=det(M), J =Aa{M){-J + QM-^Q). (13) 

Alternatively, the functions 11 and f can be derived from the topology of the corresponding 
Feynman graph G. Cutting / lines of a given connected Z-loop graph such that it becomes a 
connected tree graph T defines a chord c{T,G) as being the set of lines not belonging to this 
tree. The Feynman parameters associated with each chord define a monomial of degree I. The 
set of all such trees (or 1 -trees) is denoted by Ti. The 1 -trees T ETi define U as being the sum 
over all monomials corresponding to the chords C{T,G). Cutting one more line of a 1-tree leads 
to two disconnected trees (Ti, 72), or a 2-tree. 12 is the set of all such pairs. The corresponding 
chords define monomials of degree / + 1 . Each 2-tree of a graph corresponds to a cut defined 
by cutting the lines which connected the two now disconnected trees in the original graph. The 
square of the sum of momenta through the cut lines of one of the two disconnected trees T\ or T2 
defines a Lorentz invariant 

= { ■ (14) 

\jec(T,G) J 
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The function Jo is the sum over all such monomials times minus the corresponding invariant. 
The function f is then given by ^Fq plus an additional piece involving the internal masses nij. In 
summary, the functions U and y are obtained from the graph as follows: 



« = Ef n 



'j^C{T,G) 



{TuT2)^l2 



'j^C{TuG) 



^ = ^o + 'uY.Xjmj. (15) 

In general, 1/ is a positive semi-definite function. Its vanishing is related to the UV sub- 
divergences of the graph. Overall UV divergences, if present, will always be contained in the 
prefactor r(v — ID/ 2). In the Euclidean region, J is also a positive semi-definite function of the 
Feynman parameters xj. 

As an example we consider the graph in fig. [TJ For simplicity we assume that all internal 
propagators are massless. Then the functions U and read: 

U = X15JC23 +^15-«46+-«23-«46, 

^ = (.JCiX3.JC4+X5.JC2X6+.«lJC5^2346) (-Pl) 

+ {X(,X3X5 +X4^iX2 +X4JC6JC1235) (-P2) 

+ {X2X4X5 +X3JC1X6 +X2JC3JC1456) {-pI) ■ (16) 

Here we used the notation that Xij„,r = Xi-\-Xj -\- ... -fx,-. 

Finally let us remark, that in eq. ([7]) we restricted ourselves to scalar integrals, where the 
numerator of the integrand is independent of the loop momentum. A priori more complicated 
cases, where the loop momentum appears in the numerator might occur. However, there is a 
general reduction algorithm, which reduces these tensor integrals to scalar integrals [34,35]. The 
price we have to pay is that these scalar integrals involve higher powers of the propagators and/or 
have shifted dimensions. Therefore we considered in eq. ^ shifted dimensions and in eq. (|7]) 
arbitrary powers of the propagators. In conclusion, the integrals of the form as in eq. (|7]) are the 
most general loop integrals we have to solve. 



3 The Mellin-Barnes transformation 

In sect. [2I we saw that the Feynman parameter integrals depend on two graph polynomials U 
and J , which are homogeneous functions of the Feynman parameters. In this section we will 
continue the discussion how these integrals can be performed and exchanged against a (multiple) 
sum over residues. The case, where the two polynomials are absent is particular simple: 

nr(v,) 
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With the help of the Mellin-Barnes transformation we now reduce the general case to eq. (fTTI) . 
The Mellin-Bames transformation reads 



(Ai+A2 + ...+A„) ' = p^- — ^ / doi... [ dOn-x (18) 

— ioo — 700 

xr(-oi)...r(-o„_i)r(oi + ...+o„_i+c)A^'...A^!:i'A,r'"-"^"-'"^ 

Each contour is such that the poles of r(— o) are to the right and the poles of r(a + c) are to 
the left. This transformation can be used to convert the sum of monomials of the polynomials U 
and J into a product, such that all Feynman parameter integrals are of the form of eq. (fTTI ). As 
this transformation converts sums into products it is the "inverse" of Feynman parametrisation. 
Eq. (fTSl) is derived from the theory of Mellin transformations: Let h{x) be a function which is 
bounded by a power law for x ^ and x^°°, e.g. 

\h{x)\<Kx^''' forx^O, 

\h{x)\<K'x'' fox x^ 00. (19) 
Then the Mellin transform is defined for cq < Re o < ci by 

00 

h^{o) = Jdxh{x)x''-\ (20) 


The inverse Mellin transform is given by 

h{x) = 2^ / doh^{a)x-''. (21) 

The integration contour is parallel to the imaginary axis and cq < Re y < ci. As an example for 
the Mellin transform we consider the function 

h{x) = ^ , (22) 
with Mellin transform hy^ (a) = r(-a)r(a + c) /T{c) . For Re(-c) < Re y < we have 

1 rn^n^if),-... (23) 



(1+jc)^' Ini J r(c) 
From eq. (1231) one obtains with x = B/A the Mellin-Barnes formula 

{A+Br = — . (24) 

^ ^ 2ni J r(c) 
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Eq. (1181) is then obtained by repeated use of eq. (1241) . 

With the help of eq. (flTl) and eq. (fTSl) we may exchange the Feynman parameter integrals 
against multiple contour integrals. A single contour integral is of the form 

J ^ J_ f r(a + ai)...r(a + a^) r{-o + bi)...r{-o + bn) 

2ni J ^ r{(5+c2)...r{o+cp)r{-o+di)...r{-(5+dq)^ 

If max (Re ( —a i ),..., Re( —am) ) < min (Re (Z?i ),..., Re(Z7„) ) the contour can be chosen as a straight 
line parallel to the imaginary axis with 

max(Re(-ai),...,Re(-am)) < Rey < min(Re(Z7i), ...,Re(£'„)) , (26) 

otherwise the contour is indented, such that the residues of r(a + ai), r(a + a„,) are to the 
right of the contour, whereas the residues of r(— a + /?i), r{—(5 + b„) are to the left of the 
contour. We further set 

a = m + n — p — q, 
(3 = m — n — p + q, 

(in n p ^ \ 1 

E«j + E^7-L^i- L^/ ]-:^im + n-p-q). (27) 
7=1 7=1 J=l J=l / ^ 

Then the integral eq. (l25l) converges absolutely for a > [36] and defines an analytic function in 

|arg.x:| < min^7l,a^^. (28) 

The integral eq. (|25] ) is most conveniently evaluated with the help of the residuum theorem by 
closing the contour to the left or to the right. Therefore we need to know under which condi- 
tions the semi-circle at infinity used to close the contour gives a vanishing contribution. This is 
obviously the case for \x\ < 1 if we close the contour to the left, and for |x| > 1, if we close the 
contour to the right. The case \x\ = 1 deserves some special attention. One can show that in the 
case P = the semi-circle gives a vanishing contribution, provided 

I < -1. (29) 



To sum up all residues which lie inside the contour it is useful to know the residues of the Gamma 
function: 

(-1)" 

res (r(o + a),o= —a — n) = ; — , res {r{—o + a),o — a + n) 

nl 

In general, one obtains (multiple) sum over residues. In particular simple cases the contour 
integrals can be performed in closed form with the help of two lemmas of Barnes. Barnes first 



(-ir 

nl 

(30) 
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lemma states that 



[CD 

^ J dar{a + a)r{b + a)r{c-a)r{d-a) 



r{a + c)r{a + d)T{h + c)T{b + d) 



Ini J ^ ' ^ ' ^ ' ^ ' r(a + Z7 + c + J) 

(31) 

if none of the poles of T{a + Q)V{h + o) coincides with the ones from T{c — <5)T{d — a). Barnes 
second lemma reads 



-/■ 



loo 

r{a + o)r{b + a)r(c + a)r( J - o)r{e - o) 

r{a + b + c + d + e + (5) 

o 

r{a + d)T{b + d)T{c + d)Y{a + e)Y{b + e)T{c + e) 

r(fl+Z7+6?+e)r(a+c+i/+e)r(Z7+c+i/+e) 



(32) 



Although the Mellin-Bames transformation has been known for a long time, the method has seen 
a revival in applications in recent years [17-19, 37^9]. 



4 Shuffle algebras 

Before we continue the discussion of loop integrals, it is useful to discuss first shuffle algebras 
and generalisations thereof from an algebraic viewpoint. Consider a set of letters A. The set A is 
called the alphabet. A word is an ordered sequence of letters: 

w = hh-.-h- (33) 

The word of length zero is denoted by e. Let ^ be a field and consider the vector space of words 
over K. A shuffle algebra si on the vector space of words is defined by 

{hh-.-h) ■{lk+\---lr) = l<;{\)lo{2)---lo{r), (34) 

shuffles O 

where the sum runs over all permutations o, which preserve the relative order of 1 , 2, . . . , /: and of 
fc+ 1, ...,r. The name "shuffle algebra" is related to the analogy of shuffling cards: If a deck of 
cards is split into two parts and then shuffled, the relative order within the two individual parts is 
conserved. The empty word e is the unit in this algebra: 

e - w = w • e = w- (35) 

A recursive definition of the shuffle product is given by 

{hh-h) ■ ik+l-lr) = h [{h-h) ■ {lk+l-lr)]+lk+l [{hh-h) ■ ik+l-lr)] (36) 

It is well known fact that the shuffle algebra is actually a (non-cocommutative) Hopf algebra [50]. 
In this context let us briefly review the definitions of a coalgebra, a bialgebra and a Hopf algebra. 
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which are closely related: First note that the unit in an algebra can be viewed as a map from K 
to A and that the multiplication can be viewed as a map from the tensor product A®AX.oA (e.g. 
one takes two elements from A, multiplies them and gets one element out). 

A coalgebra has instead of multiplication and unit the dual structures: a comultiplication A 
and a counit e. The counit is a map from A to K, whereas comultiplication is a map from AioA® 
A. Note that comultiplication and counit go in the reverse direction compared to multiplication 
and unit. We will always assume that the comultiplication is coassociative. The general form of 
the coproduct is 

A(a) = (37) 

where af^ denotes an element of A appearing in the first slot of A® A and af^ correspond- 
ingly denotes an element of A appearing in the second slot. Sweedler's notation [51] consists in 
dropping the dummy index / and the summation symbol: 

A(a) = a(^)®a(^) (38) 

The sum is implicitly understood. This is similar to Einstein's summation convention, except 
that the dummy summation index i is also dropped. The superscripts and indicate that a 
sum is involved. 

A bialgebra is an algebra and a coalgebra at the same time, such that the two structures are 
compatible with each other. Using Sweedler's notation, the compatibility between the multipli- 
cation and comultiplication is expressed as 

A(a-Z7) = (a(i)-Z7(i))® (a(2).^(2)j_ (39^ 

A Hopf algebra is a bialgebra with an additional map from A to A, called the antipode 5 , 
which fulfils 

a(^)-5 =5 (a(^))-a(2)=0 foxai^e. (40) 

With this background at hand we can now state the coproduct, the counit and the antipode 
for the shuffle algebra: The counit e is given by: 

e{e) = \, e{kh...ln)=0. (41) 

The coproduct A is given by: 

k 

Aikh-k) = £(/;■+!... 4) (42) 

7=0 

The antipode S is given by: 

S{hl2...lk) = {-If hh-x-hh- (43) 
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Figure 2: Sketch of the proof for the shuffle product of two iterated integrals. The integral over 
the square is replaced by two integrals over the upper and lower triangle. 



The shuffle algebra is generated by the Lyndon words. If one introduces a lexicographic ordering 
on the letters of the alphabet A, a Lyndon word is defined by the property 

w <v (44) 

for any sub-words u and v such that w = uv. 

An important example for a shuffle algebra are iterated integrals. Let [a,b\ be a segment of 
the real line and fi, fi, ... functions on this interval. Let us define the following iterated integrals: 

b t\ h-i 

I{fij2,-Jk;a,b) = J Mti)dh J fi{ti)dt2... J fk{tk)dtk (45) 

a a a 

For fixed a and b we have a shuffle algebra: 

Hfhf2,-,fk;a,b)-I{fk+i,...,fr;a,b) = /(/(,(i),/(^(2), a, Z?), (46) 

shuffles O 

where the sum runs over all permutations a, which preserve the relative order of 1,2, and 
of k+ I, ...,r. The proof is sketched in fig. [2l The two outermost integrations are recursively 
replaced by integrations over the upper and lower triangle. 

We now consider generalisations of shuffle algebras. Assume that for the set of letters we 
have an additional operation 

(.,.) : A®A^A, 

h®li^{hM). (47) 

which is commutative and associative. Then we can define a new product of words recursively 
through 

{hh-h) * ik+i-h) = h [{li--h) * {lk+\--U)]+lk+\ [{hh-h) * {k+i-h)] 

+ihJk+i)[il2...lk)Hlk+i...lr)] (48) 

This product is a generalisation of the shuffle product and differs from the recursive definition 
of the shuffle product in eq. (|36l ) through the extra term in the last line. This modified product 
is known under the names quasi-shuffle product [52], mixable shuffle product [53] or stuffle 
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Figure 3: Sketch of the proof for the quasi-shuffle product of nested sums. The sum over the 
square is replaced by the sum over the three regions on the r.h.s. 



product [54]. Quasi-shuffle algebras are Hopf algebras. Comultiplication and counit are defined 
as for the shuffle algebras. The counit e is given by: 

e{e) = h eihl2...ln)=0. (49) 

The coproduct A is given by: 

k 

A{hl2...lk) = j;(/;+1...4)®(/l. ../;). (50) 

The antipode S is recursively defined through 

k-l 

sikh-k) = -hl2...lk-Y.^{h+i-h)*{h...h). (51) 

7=1 

An example for a quasi-shuffle algebra are nested sums. Let and be integers with < 
and let f\, fi, ... be functions defined on the integers. We consider the following nested sums: 

S{fuf2,-Jk;na,nt) = J2 L Mi2)-- fkiik) (52) 

For fixed na and nt, we have a quasi-shuffle algebra: 

S{f\j2,---Jk\na,nb)*S{fk^i,...Jr\na,ni,) = 

ni, 

S{f2-,-Jk\na,il - l)*S{fk+U...Jr;naJl - 1) 

i\=na 
lib 

+ J2 fkUl)S{flj2,---Jk\naJl-l)*S{fk+2-,---Jr\naJl-l) 

i\=na 

+ fi{i)fk{i) S{f2-,---Jk,na-,i-\)*S{fk+2-:---Jr;naJ-\) (53) 

Note that the product of two letters corresponds to the point- wise product of the two functions: 

{fufj){n) = Mn)fj{n). (54) 

The proof that nested sums obey the quasi-shuffle algebra is sketched in Fig. |3l The outermost 
sums of the nested sums on the l.h.s of (l53l) are split into the three regions indicated in Fig. |3l 
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5 Multiple polylogarithms 



In the previous section we have seen that iterated integrals form a shuffle algebra, while nested 
sums form a quasi-shuffle algebra. In this context multiple polylogarithms form an interesting 
class of functions. They have a representation as iterated integrals as well as nested sums. There- 
fore multiple polylogarithms form a shuffle algebra as well as a quasi-shuffle algebra. The two 
algebra structures are independent. Let us start with the representation as nested sums. The 
multiple polylogarithms are defined by 



^1 

/i>/2>...>/;t>0 1 ^ 



The multiple polylogarithms are generalisations of the classical polylogarithms Li„(jc) [55], 
whose most prominent examples are 



Lil W = E — = -ln(l -x), U2{x) = £ r^, (56) 



(1=1 ^1 (,=1 h 



as well as Nielsen's generalised polylogarithms [56] 



Sn,p{x) = Li„+i,i,...^i(x,l^^^, (57) 

p-l 

and the harmonic polylogarithms [57] 

Hm\,...,mk{x) = himi,...,mii{x, i, 1 )- (58) 

k-l 

Multiple polylogarithms have been studied extensively in the literature by physicists [57-70] and 
mathematicians [54,71-81]. 

In addition, multiple polylogarithms have an integral representation. To discuss the integral 
representation it is convenient to introduce for ZkT^O the following functions 

y ti 

n( \ fdtifdh f dtk 

G{zi,...,Zk;y) = /; / ; • (59) 

J h—z\ J ti — zi J tk — Zk 



In this definition one variable is redundant due to the following scaling relation: 

G{zh...,Zk;y) = G{xzh...,xzk;xy) (60) 

If one further defines 

8{z;y) = (61) 
y-z 
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then one has 



^G{zi,...,Zk;y) = g{zi;y)G{z2,---,Zk;y) (62) 



and 



G(zi,z2,---,Zyt;j) = J dt g{zi;t)G{z2,---,Zk;t). 



(63) 



One can slightly enlarge the set and define G(0, ...,0; j) with k zeros for zi to Zk to be 

G(0,...,0;y) = l^i^nyf. (64) 

This permits us to allow trailing zeros in the sequence (zi , Zjt) by defining the function G with 
trailing zeros via (l63l) and (|64l) . To relate the multiple poly logarithms to the functions G it is 
convenient to introduce the following short-hand notation: 

Gmi,-,mk{zi,...,zk;y) = G{0^^^^,zi,...,zk-i,0^^^,zk;y) (65) 
Here, all for j = 1, are assumed to be non-zero. One then finds 

Limi,...,m^(-^l, — ( — l)^^^!,...,^^. ( — , Uj- (66) 



Xl X\X2 Xl...Xk 



The inverse formula reads 



Gmx,...,mk{Zi,...,Zk\y) = (-l)'^Limi,...,mJ ) • (67) 

\Zl Z2 Zk J 

Eq. (|66l ) together with (|65] ) and (|59l ) defines an integral representation for the multiple polyloga- 
rithms. To make this more explicit I first introduce some notation for iterated integrals 



A A t„ t2 

J t — Un t — ai J tn—anJ tn-\—an-\ J 





and the short hand notation: 



(68) 



A , A 

dt \ dt f dt dt dt 

ol = —0... — . (69) 

t — a J J ^ t ^ t — a 

m times 
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The integral representation for himi,...,mi,{xi, ■■■,Xk) reads then 



^ dt dt 



t-b 



dt dt fdt \""^-^ dt 



o ^o...o -o -, (70) 

t J t-b2 \t J t-bk 

where the bfs are related to the xfs 

bj = i . (71) 

XlX2...Xj 

Up to now we treated multiple polylogarithms from an algebraic point of view. Equally 
important are the analytical properties, which are needed for an efficient numerical evaluation. 
As an example I first discuss the numerical evaluation of the dilogarithm [82]: 



X 



, , ln(l-0 f x« 

U2ix) = -jdt^^=Z-, 



(72) 







The power series expansion can be evaluated numerically, provided |x| < 1. Using the functional 
equations 



U2{x) = -Li2Q)-|-^(ln(-x))2, 



U2{x) = -Li2(l-A;) + — -ln(jc)ln(l-A;). (73) 

o 

any argument of the dilogarithm can be mapped into the region \x\ < 1 and — 1 < Re(x) < 1 /2. 
The numerical computation can be accelerated by using an expansion in [— ln(l — x)] and the 
Bernoulli numbers Bf. 

CO n 

Li2W = Itt-W^-MI-^))'^'- (74) 

The generalisation to multiple polylogarithms proceeds along the same lines [67]: Using the 
integral representation 

Gmu...,mk {ZI,Z2, ■■■,Zk;y) = (75) 



t — Z\\t J t—Z2 \t I t—Zk 



15 



one transforms all arguments into a region, where one has a converging power series expansion: 

°° °° 1 / -1, \ -/i 



°° °° ^ 3^ A ' 

Gmu...,mk{zh...:Zk;y) = ••H T~Z — I wi ( — 



1 /jV^ 1 /yV^ 



The multiple poly logarithms satisfy the Holder convolution [54]. For zi 7^ 1 and Zw 7^ this 
identity reads 

G(zi,...,Zvv;l) = (77) 
(-l)^G 1^1 -zy, 1 1 -zi; 1 - -^^ G l^z^+i, ...,Zw; . 

The Holder convolution can be used to accelerate the convergence for the series representation 
of the multiple poly logarithms. 



6 Laurent expansion of Feynman integrals 

Let us return to the question on how to compute Feynman integrals. In section [3] we saw how to 
obtain from the Mellin-Bames transformation (multiple) sums by closing the integration contours 
and summing up the residues. As a simple example let us consider that the sum of residues is 
equal to 

r(/+l)r(/ + a3 + ?3e) 

Here ai, ^2 and a-i are assumed to be integers. Up to prefactors the expression in eq. (1781) is a 
hyper-geometric function 2F\ . We are interested in the Laurent expansion of this expression in 
the small parameter £. The basic formula for the expansion of Gamma functions reads 

r(n + e) =r(l+£)r(n) [\+zZi{n-\)+£^Zn{n-l) 

WZni{n - 1) + ... + £"-iZii...i(n - 1)] , (79) 

where Zmi,...,mi,i>T^) are Euler-Zagier sums defined by 

Zmu...,m,{n) = £ Ai---A^- (80) 

n>ii>i2>...>ik>0^ ^ 

This motivates the following definition of a special form of nested sums, called Z-sums: 

Z(n;mi,...,mk;xi,...,Xk) = (81) 
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k is called the depth of the Z-sum and w = mi + ... + m^;- is called the weight. If the sums go to 
infinity (n = oo) the Z-sums are multiple poly logarithms: 

Z{oo;mi,...,mk,xi,...,Xk) = Li„,,...,„^(xi, ...,JC;t)- (82) 
For xi = ... = Xk = I the definition reduces to the Euler-Zagier sums [83, 84]: 

Z(n;mi,...,m;t;l, •••,!) = Z„Ji,...,mi («)■ (83) 

For n = oo and xi = ... = Xk = I the sum is a multiple ^- value [54]: 

Z(oo;mi,...,mfc; !,...,!) = (84) 

The usefulness of the Z-sums lies in the fact, that they interpolate between multiple polyloga- 
rithms and Euler-Zagier sums. The Z-sums form a quasi- shuffle algebra. 

Using r{x+ 1) = xr{x), partial fractioning and an adjustment of the summation index one 
can transform eq. (fTSl) into terms of the form 

y ni + tie)r{i + t2e) X' 



where m is an integer. Now using eq. ([791) one obtains 

r^^. , (l+e^iZi(/-l) + -)(l+e^2Zi(/-l)-f...) ^' 

"^^'+'1^1 (l+e^3Zi(/-l) + ...) ^''^ 

Inverting the power series in the denominator and truncating in £ one obtains in each order in £ 
terms of the form 

^Z„,,. ..,„,(/- 1)Z„/ 1)Z„// 1) (87) 

!=1 ^ 

Using the quasi-shuffle product for Z-sums the three Euler-Zagier sums can be reduced to single 
Euler-Zagier sums and one finally arrives at terms of the form 

oo / 

£^Z„,..,„,(/-1), (88) 

which are special cases of multiple polylogarithms, called harmonic polylogarithmsi/m^,„j^...^,„^(x). 
This completes the algorithm for the expansion in £ for sums of the form as in eq. (TTSl) . 

The Hopf algebra of Z-sums has additional structures if we allow expressions of the form 

-^Z{n\mi,...,mk\xi,...,Xk), (89) 



e.g. Z-sums multiplied by a letter. Then the following convolution product 

n — 1 / n — 7 

r y 



I ];;^20--l;-)T^Z(n-/-l;...) (90) 
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can again be expressed in terms of expressions of the form (|89| ). In addition there is a conjugation, 
e.g. sums of the form 



i J(-i)'i^20;-) (91) 



can also be reduced to terms of the form (1891) . The name conjugation stems from the following 
fact: To any function f{n) of an integer variable n one can define a conjugated function C o f{n) 
as the following sum 

Cof{n) = i;(^)(-l)7(0- (92) 

Then conjugation satisfies the following two properties: 

Col = 1, 

CoCofin) = f{n). (93) 
Finally there is the combination of conjugation and convolution, e.g. sums of the form 



-E(,)(-»'^Z(,;...)^Z(„-/;...) (94) 



can also be reduced to terms of the form (1891) . These properties can be used to expand more 
complicated transcendental functions like 



oo oo 



or 



^ ^ Tji + ai) Tji + ak) Tjj + h) Tjj + bi)r{i + j + ci) Tji + j + c^) , j 

ro;tor(/+a;)--T(/+4)r(j>z.;)--T(j>z.;)r(/+j>c;)--T(/+j>o^^ ^ ^ 



i+j\r{i+ai) r{i+ak)TU+bi) r{j +bi)r{i+ j +ci) r{i+j+c,n)_ ^ j 



hh^ j ;r(/+a;)-r(/+4) r{j+b[y-nj+b',)r{i+j+c[yT{i+j+c'j''^' 

(96) 

Examples for functions of this type are the first and second Appell function Fi and F2. Note 
that in these examples there are always as many Gamma functions in the numerator as in the 
denominator. We assume that all a„, a'„, b„, b'„, Cn and c', are of the form "integer + const ■ £". 

The first type can be generalised to the form "rational number + const ■ £", if the Gamma 
functions always occur in ratios of the form 

T{n + a-^+bz) 

1 (97) 

r(n + c-f+Je)' ^ ^ 
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where the same rational number p/q occurs in the numerator and in the denominator [66]. In 
this case we have to replace eq. (1791 ) by 



r^n + l-^ + £ 



r(i-£+e)r(n + i-g 



) 



(98) 





which introduces the q-th roots of unity 



exp 




(99) 



In summary these techniques allow a systematic procedure for the computation of Feynman 
integrals, if certain conditions are met. These conditions require that factors of Gamma functions 
are balanced like in eq. (|95l) or eq. (|96l) [63, 66]. The algebraic properties of nested sums and 
iterated integrals discussed here are well-suited for an implementation into a computer algebra 
system and several packages for these manipulations exist [58, 85-88]. 

7 Conclusions 

In this article I discussed the mathematical structures underlying the computation of Feynman 
loop integrals. One encounters iterated structures as nested sums or iterated integrals, which form 
a Hopf algebra with a shuffle or quasi-shuffle product. Of particular importance are multiple 
poly logarithms. The algebraic properties of these functions are very rich: They form at the same 
time a shuffle algebra as well as a quasi-shuffle algebra. Based on these algebraic structures I 
discussed algorithms which evaluate Feynman integrals to multiple poly logarithms. 
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